Search results for "Mathematics - Numerical Analysis"
showing 10 items of 87 documents
Stochastic Galerkin method for cloud simulation
2018
AbstractWe develop a stochastic Galerkin method for a coupled Navier-Stokes-cloud system that models dynamics of warm clouds. Our goal is to explicitly describe the evolution of uncertainties that arise due to unknown input data, such as model parameters and initial or boundary conditions. The developed stochastic Galerkin method combines the space-time approximation obtained by a suitable finite volume method with a spectral-type approximation based on the generalized polynomial chaos expansion in the stochastic space. The resulting numerical scheme yields a second-order accurate approximation in both space and time and exponential convergence in the stochastic space. Our numerical results…
Better numerical approximation by Durrmeyer type operators
2018
The main object of this paper is to construct new Durrmeyer type operators which have better features than the classical one. Some results concerning the rate of convergence and asymptotic formulas of the new operator are given. Finally, the theoretical results are analyzed by numerical examples.
Estimates for the differences of positive linear operators and their derivatives
2019
The present paper deals with the estimate of the differences of certain positive linear operators and their derivatives. Oxur approach involves operators defined on bounded intervals, as Bernstein operators, Kantorovich operators, genuine Bernstein-Durrmeyer operators, and Durrmeyer operators with Jacobi weights. The estimates in quantitative form are given in terms of the first modulus of continuity. In order to analyze the theoretical results in the last section, we consider some numerical examples.
Better approximation of functions by genuine Bernstein-Durrmeyer type operators
2018
The main object of this paper is to construct a new genuine Bernstein-Durrmeyer type operators which have better features than the classical one. Some direct estimates for the modified genuine Bernstein-Durrmeyer operator by means of the first and second modulus of continuity are given. An asymptotic formula for the new operator is proved. Finally, some numerical examples with illustrative graphics have been added to validate the theoretical results and also compare the rate of convergence.
Rough McKean-Vlasov dynamics for robust ensemble Kalman filtering
2021
Motivated by the challenge of incorporating data into misspecified and multiscale dynamical models, we study a McKean-Vlasov equation that contains the data stream as a common driving rough path. This setting allows us to prove well-posedness as well as continuity with respect to the driver in an appropriate rough-path topology. The latter property is key in our subsequent development of a robust data assimilation methodology: We establish propagation of chaos for the associated interacting particle system, which in turn is suggestive of a numerical scheme that can be viewed as an extension of the ensemble Kalman filter to a rough-path framework. Finally, we discuss a data-driven method bas…
An elementary formula for computing shape derivatives of EFIE system matrix
2012
We derive analytical shape derivative formulas of the system matrix representing electric field integral equation discretized with Raviart-Thomas basis functions. The arising integrals are easy to compute with similar methods as the entries of the original system matrix. The results are compared to derivatives computed with automatic differentiation technique and finite differences, and are found to be in excellent agreement.
Guaranteed error control bounds for the stabilised space-time IgA approximations to parabolic problems
2017
The paper is concerned with space-time IgA approximations of parabolic initial-boundary value problems. We deduce guaranteed and fully computable error bounds adapted to special features of IgA approximations and investigate their applicability. The derivation method is based on the analysis of respective integral identities and purely functional arguments. Therefore, the estimates do not contain mesh-dependent constants and are valid for any approximation from the admissible (energy) class. In particular, they provide computable error bounds for norms associated with stabilised space-time IgA approximations as well as imply efficient error indicators enhancing the performance of fully adap…
Existence of dynamical low-rank approximations to parabolic problems
2021
The existence and uniqueness of weak solutions to dynamical low-rank evolution problems for parabolic partial differential equations in two spatial dimensions is shown, covering also non-diagonal diffusion in the elliptic part. The proof is based on a variational time-stepping scheme on the low-rank manifold. Moreover, this scheme is shown to be closely related to practical methods for computing such low-rank evolutions.
A degenerating convection-diffusion system modelling froth flotation with drainage
2022
Abstract Froth flotation is a common unit operation used in mineral processing. It serves to separate valuable mineral particles from worthless gangue particles in finely ground ores. The valuable mineral particles are hydrophobic and attach to bubbles of air injected into the pulp. This creates bubble-particle aggregates that rise to the top of the flotation column where they accumulate to a froth or foam layer that is removed through a launder for further processing. At the same time, the hydrophilic gangue particles settle and are removed continuously. The drainage of liquid due to capillarity is essential for the formation of a stable froth layer. This effect is included into a previous…
A generalized Newton iteration for computing the solution of the inverse Henderson problem
2020
We develop a generalized Newton scheme IHNC for the construction of effective pair potentials for systems of interacting point-like particles.The construction is made in such a way that the distribution of the particles matches a given radial distribution function. The IHNC iteration uses the hypernetted-chain integral equation for an approximate evaluation of the inverse of the Jacobian of the forward operator. In contrast to the full Newton method realized in the Inverse Monte Carlo (IMC) scheme, the IHNC algorithm requires only a single molecular dynamics computation of the radial distribution function per iteration step, and no further expensive cross-correlations. Numerical experiments…